function h = horizontal_d(i)
    h = i(:, 2:end) - i(:, 1:end-1);
end